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Abstract 

We include p-spin interactions in a spherical version of a soluble mean-field spin- 
glass model proposed by van Hemmen. Due to the simplicity of the solutions, which 
do not require the use of the replica trick, we are able to carry out a detailed 
investigation of a number of special situations. For p > 3, there appear first-order 
transitions between the paramagnetic and the ordered phases. In the presence of 
additional ferromagnetic interactions, we show that there is no stable mixed phase, 
with both ferromagnetic and spin-glass properties. 



1 Introduction 



The attempts to account for the physical properties of disordered systems have 
been a continuous source of models and problems in statistical mechanics. It is 
known that a mean-field p-spin-glass model, with p > 3, displays a first-order 
transition to a glassy phase, and non-trivial dynamical behavior. Although 
models with p-spin interactions may be far from representing real glasses, 
it has been noted that the dynamical mean-field equations do resemble the 
corresponding mode-coupling expressions for glassy systems [1]. 

In some recent publications, a ferromagnetic term, as well as additional multi- 
spin interactions, have been added to the simple spin-glass p-spin model [2,3,4]. 
These new systems display continuous and first-order transitions, and a glassy 
ferromagnetic region with pronounced reentrant borders. However, the search 
for the thermodynamic solutions of these mean-field models requires the use 
of the subtleties of the replica method. The replica-symmetric solutions have 
to be supplemented by a one-step replica-symmetry breaking scheme in order 
to characterize the glassy ferromagnetic phase. 

In this paper, we decided to investigate the thermodynamic behavior of highly 
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simplified versions of disordered multi-spin mean-field models. The idea con- 
sists in the inclusion of p-spin interactions, with p > 3, in a slightly more 
general version of a spin-glass model proposed by van Hemmen about twenty 
years ago [5]. This simple van Hemmen model, which does not require the use 
of the replica method, is obviously unable to reproduce the rich (ultramet- 
ric) structure of the spin-glass phase [6]. However, it does give a reasonable 
account of the main features of the phase diagrams, including the location 
of phase regions and transition lines. In this paper, we further simplify the 
calculations by assuming spherical instead of Ising spin variables. 

In Section 2, we review the van Hemmen spin-glass model with spherical spin 
variables. In Section 3, we introduce multi-spin interactions in the van Hem- 
men model, and make a number of contacts with recent work. Section 4 is 
devoted to some conclusions. 



2 The van Hemmen model 



The van Hemmen spin-glass model, with the inclusion of a ferromagnetic term, 
is given by the Hamiltonian 

J Q N 

H = -— a i a j~ Yl JijViVj - Hj2 a i, (!) 

l<i<j<N l<i<j<N i=l 

where J > 0, 

J ll = ±J^(f + ^], (2) 

J > 0, and H is an external field. The set of independent, identically dis- 
tributed random variables {^f^j, for % = 1,2, ...,N and a — 1,2, associated 
with a double-delta probability distribution, 

p(d" , ) = ^(^ , - 1 ) + ^(^' + 1 ). (») 

is supposed to mimic the presence of disorder and competition in real spin 
glasses. In the spherical version of this model, the spin variables are real num- 
bers, — oo < <jj < oo, for all i = 1,2, ...,N, with the spherical constraint, 

N 

T, a i= N - ( 4 ) 



Using more convenient variables, 

Y n ^ N 

m = T7 £ <*i, m a = T7 E ( 5 ) 
iV i=l iV i=l 



2 



for a — 1,2, and discarding irrelevant terms in the thermodynamic limit, we 
can rewrite the Hamiltonian, given by Eq. (1), in the form 



1 9 1 
H = --J Nm 2 - -JNm 1 m 2 - HNm. 



(6) 



Given a configuration of the random variables, the partition function may be 
written as 



Z = Tr exp 



1 1 

-pJ Nm 2 + -pJNm 1 m 2 + (3 HNm 



(7) 



where (5 is the inverse of temperature, and the trace should take into account 
the spherical constraint. 

We now rewrite the partition function in the form 



/+oo r+oo r+oo 

dm / dm\ / dm 2 il (m,mi,m 2 ) 
-oo J — oo J — oo 

x exp J Q Nm 2 + ^(3JNm 1 m 2 + pHNmj , 



where 



(8) 



"= mi (£» 5 (j> 2 - ") 5 (— ^g>) 



(9) 



In the thermodynamic limit, we use an integral representation for the delta 
functions and invoke the law of large numbers in order to obtain the asymptotic 
result [7] 



Q ~ exp 1 + In (2ir) + In (l — m 2 — m\ — m\?j 



which leads to the partition function 



/+oo r+oo r+oo 

dm j dmi I dm 2 exp [— (3Nf (to, TOi, m 2 )] , 
-oo J —oo J —oo 



(10) 



where 



m 2 — ml — m 



/ = -- ^J m 2 + Jmim 2 ) -Em-— 1 + In (2tt) + In (l 

(12) 

The thermodynamic solutions come from the stable minima of this free-energy 
functional. 
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Fig. 1. Phase diagram of the spherical van Hemmen model. Solid lines correspond 
to zero external field, while dashed lines are obtained for Hr/J = 1/4. Thin (thick) 
lines indicate second (first) order transitions. The labels correspond to the param- 
agnetic (PARA), ferromagnetic (FERRO), and to the "spin-glass" (SG) phases. 

In zero external field, H = 0, it is easy to draw the phase diagram of Figure 
1, in terms of T = (/3J) _1 and r = Jo/ J. Besides the disordered paramagnetic 
phase {m — m\ — m 2 = 0), there is a "spin-glass" (m = 0, and mi,m 2 ^ 0) 
and a ferromagnetic phase (m ^ 0, and m\ = m 2 = 0). In the ferromag- 
netic region, m 2 = 1 — T/r; in the spin-glass region, ml = m\ = 1/2 — T. 
The paramagnetic borders are lines of continuous phase transitions. Mixed 
phases (with m ^ and m 1 ,m 2 ^ 0) are restricted to a line of coexis- 
tence at r = J / J = 1/2 (which also corresponds to the limits of stability 
of the ferromagnetic and spin-glass solutions). Similar topological features are 
also present in the phase diagram of the much more elaborate Sherrington- 
Kirkpatrick (SK) mean-field model of a spin glass. 

In the presence of a random field, the Hamiltonian of the van Hemmen model 
is written as 

Jo N 
H = -— Yl JijViVj -J^Hicn, (13) 

l<i<j<N l<i<j<N i=l 

where {Hi}, for i = 1,2, ...,7V, is a set of quenched independent, identically 
distributed random variables, given by the probability distribution 

Ph m = ^5 - H R ) + ^5 (Hi + H R ) . (14) 
Introducing the new variable 



it is straightforward to obtain the free-energy functional 



/ = ~9 ( J o m2 + ^1^2) -HrV-t^ I 1 + ln ( 27 + ln (l 



222 
m — m l — m 2 



(16) 

In the T — r phase diagram, there is a depression of the paramagnetic lines, 
which still meet at r = 1/2 (see Figure 1). For small values of Hr = Hr/ J « 
1, we have the following asymptotic forms of the paramagnetic critical lines: (i) 
T = 1/2— 2h 2 R: at the spin-glass border; (ii) T = r—h 2 R /r, at the paramagnetic- 
ferromagnetic border. There is no tricritical point with spherical spin variables. 
The same qualitative depression of the paramagnetic borders is also present 
in the phase diagram of an SK model in a random field [8]. 



3 Multi-spin interactions 



The inclusion of p-spin interactions, with p > 2, leads to a natural generaliza- 
tion of the van Hemmen model. Let us write the spin Hamiltonian 



H = 



with 



Ji 



l<il<l 2 <...<%,<iV 



Jr) 



Ji\...i p ®i\ ■ "Pip 1 



'v \ " A<xi)A°u) A a p) 

jj\]\[p-l >u ■■■>>ip 1 



H... ip 



(17) 



(18) 



where the sum is over the permutations of the p patterns, a±, a p . In analogy 
with the p = 2 case, we have enlarged the set of random variables, jd"''}) 
which are still given by the probability distribution of Eq. (3). With this 
choice of patterns, the p-spin Hamiltonian can be written as 



H = — -J v Nm\...m p , 
p\ 

from which we obtain the free-energy functional 



(19) 



f 1 T 1 



1 + In (2tt) + In 1 - ^ 



(20) 



For p > 3, it is easy to show that there is a first-order transition between a 
paramagnetic disordered phase (mi = m 2 = ... = 0) and an ordered "spin- 



glass" phase (ml 



... 7^ 0). The existence of a first-order transition in 



these p-spin models, which leads to spinodal lines and may be responsible for 
peculiar dynamical phenomena, is usually taken as an important contact with 
the behavior of real glasses. 
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These calculations can be easily extended to a much larger class of models. Let 
us consider blocks of p-spin and r-spin interactions, as well as ferromagnetic 
n-spin interactions and a random external field. A sufficiently general spin 
Hamiltonian may be written as 



H = - 



E 



j, 



ii...ipvi± ■ 



■CTi„ - 



E 



j, 



l\...lr u l\ ■ 



■Oi — 



l<i 1 <i 2 <...<i p <N 



l<i 1 <i 2 <-<i r <N 



with 



and 



Jo 



n \JS[n-l 



N 



E a * 

i=l 



N 



- H i a h 



i=l 



j Jp 

J ^ p \ N p- 



Et(ai)t(aa) C^p) 



(ai...a p ) 



Ji 



^l...^ r 



r!A rr-i 



EA a P+l) A a P+2) A a P+r) 



(21) 



(22) 



(23) 



(a p +i...a p + r ) 

where the sums are over all permutations of a±...a p and a p+ i...a p+r , and the 
independent random variables are still given by the probability distributions 
of equations (3) and (14). 



Using the definitions of m, m a , and q, given by Eqs. (5) and (15), for a = 
1, ...,p+r, it is easy to rearrange the Hamiltonian in the more convenient form 



11 1 

H = — -J p Nmi...m p -NJ r m p+ i...m p+r -J Nm n — H R Nq. (24) 

In the thermodynamic limit, assuming finite values of the parameters p, r, 
and n, we can write the free-energy functional 



/ = -—J p m 1 ...m p - — J r m p+1 ...m p+r - —J Q m n - H R q 
p\ r! n\ 



1 

2~P 



p+r 

1 + In (2vr) + In j 1 - m 2 - q 2 - ^ m 2 a 

a = l 



(25) 



We now consider some particular situations: 



(A) A model with (uniform) ferromagnetic interactions J02 and J04 involving 
blocks of n = 2 and n = 4 spins, in the presence of a random field, leads to 
the free-energy functional 



1 + In (2tt) + In (l - m 2 - g 2 )] . (26) 
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1.6 




Fig. 2. Phase diagram of the mean-field spherical model with ferromagnetic J02 and 
J04 couplings in zero external field. Thin (thick) lines indicate second (first) order 
transitions, while P± r marks the tricritical point. 



T 




r 12.0 



Fig. 3. Phase diagram of the mean-field spherical model with ferromagnetic J02 
and J04 couplings in a random field of strength Hr. The surface corresponds to a 
paramagnetic boundary, and the thick curve is a line of tricritical points, separat- 
ing a region of first-order (large r) from a region of second-order (small r) phase 
transitions. 

In zero random field, Hr = 0, the phase diagram, in terms of temperature, 
T = (/3Jo2)~ ; versus the ratio of ferromagnetic uniform interactions, r = 
J04/ J02, displays a line of second-order transitions that turns into a first-order 
boundary at a tricritical point, J04/J02 = 6 (see Figure 2). The existence of 
first and second-order transitions at the mean-field level gives an additional 
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motivation to pursue the investigation of the spin-glass versions of this model. 
The random field introduces just a shift of the tricritical point. The global 
phase diagram, in terms of T, J 04 / J 02 , and the strength of the random field, 
H R /J 02 , as shown in Figure 3, displays a line of tricritical points. 

(B) If we assume r = and n = 2, in zero field, it is possible to make contact 
with some recent calculations for a spherical version of a p-spin Sherrington- 
Kirkpatrick (SK) spin-glass model with the addition of a simple ferromagnetic 
term. According to these calculations [2,3], for p > 3, the phase diagram of the 
spherical p-spin SK model displays a glassy ferromagnetic phase, besides the 
usual paramagnetic, spin glass, and pure ferromagnetic phases. Also, there is 
a pronounced reentrance of the border between the two ferromagnetic phases. 

In the context of the van Hemmen models, we write the free-energy functional 



* 1 7 1 7 2 1 

f =- — J p m 1 ...m p - -J m - — 



1 + In (27r) + In 1 - m 2 - ^ ml 



(27) 

We now set p = 3, which gives typical results for p > n. The minimization 
of / leads to the phase diagram shown in Figure 4, in terms of temperature, 
T = 1/ (PJ3), and the ratio of interactions, r = J / J3. The paramagnetic phase 
(m — mi — m 2 = m 3 = 0) is stable for T > r. The ferromagnetic phase (m 2 = 
1 — T/r, m 1 — m 2 = m 3 = 0) is stable for T < r. The critical line between 
the paramagnetic and the ferromagnetic phases ends at a critical endpoint (at 
T = r = 0.03278...). The "spin-glass phase" (m = 0, m\ = m\ = m\ 7^ 0) 
is limited by a first-order boundary (at T = 0.03278..., for small values of r; 
and at r = l/9\/3 for T = 0). The extremization of the free energy / also 
leads to a "mixed ferromagnetic solution", m^0, with m\ = m\ = m\ 7^ 0, 
which might be the analogue of the glassy ferromagnetic phase of the p-spin 
SK models. However, it is easy to show that the free energy of this solution is 
always larger than the free energy of the uniform ferromagnetic phase. For this 
p = 3 example, at any temperature, / (ferro) — / (mixed ferro) = — 18r 3 < 0. 
Indeed, it is easy to obtain analogous results for p > 3. 



4 Conclusions 



We have introduced multi-spin interactions in the van Hemmen spin-glass 
model with spherical spin variables. Using standard techniques, we write ana- 
lytic expressions for a self-averaged free energy of a model Hamiltonian includ- 
ing blocks of disordered and uniform interactions, in the presence of random 
fields. From this free energy, it is easy to draw some phase diagrams, in order 
to compare with results for the corresponding SK models. 
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Fig. 4. Phase diagram of the spherical van Hemmen model with ferromagnetic pair 
interactions and random p = 3 couplings. Thin (thick) lines indicate second (first) 
order transitions. 

In particular, for p = 3 and two-spin uniform ferromagnetic interactions, 
we obtain a phase diagram with first-order transitions between spin-glass 
and either paramagnetic or ferromagnetic phases, and a second-order ferro- 
paramagnetic transition line, which ends at a critical endpoint. In contrast 
to the calculations for the analogous SK model, we show that there is no 
possibility of appearance of a ferromagnetic phase of spin-glass character. 
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